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The coupled rjN , ttN, ^N system is described by a K-matrix method. The parameters in this 
model are adjusted to get an optimal fit to ttN -^ tvN, nN -^ rjN, jN -^ tvN and 7A'' -^ rjN 
data in an energy range of about lOOMeV or so each side of the 77 threshold. The outcome is the 
appearance of two solutions one which has an »7-nucleon scattering length (a) of about 1.0 fm and 
a second with a ~ 0.2fm. However, the second solution has an unconventional non-Lorentian form 
f^ ' for the T-matrix in the region of the 1535(20)MeV and 1650(30)MeV S-wave ttN resonances. 

o 

CN ■ PACS numbers: 13.75.-n, 25.80.-e, 25.40.V 



o 



The value of the 77-nuclcon scattering length (a) is still uncertain, but with everyone agreeing that it is indeed 
attractive i.e. a > 0. In the literature estimates can be found ranging from Re a < 0.3fm ||l[ upto about l.OiO.lfm 
- a selection being given in Table I. 

The main interest in a lies in the fact that, if a is sufficiently attractive, then ry-nuclear quasi-bound states may 



be possible. These were first suggested about 15 years ago |12 13]. Since then many articles have appeared on this 

I , subject studying different reactions in which such quasi-bound states could manifest themselves e.g. in the 77-deuteron 

J> ' system [Q, the pd -^ ^Herj reaction ||l^,|l6| and in the pp -^ pprj cross section ||l^. However, the first experimental 

C*n ' attempt to discover restates in heavier nuclei gave a negative conclusion |Q. Another experiment is now proposed 

10 to check this result [^,Q. Unfortunately, in the absence of 77-beams, these same reactions are also the only source of 

^^ ■ experimental information about the rjN scattering length. It is, therefore, important that in any discussion as many 

reactions as possible are treated simultaneously. Otherwise success with one reaction may be completely nullified by 

failure with another. 



o 

(<— ») , With this in mind, in Ref. |^1[], the present authors carried out a simultaneous iiT-matrix fit to the nN -^ rjN 

■^ ■ cross sections reviewed by Ncfkcns p3| and the 77? -^ rjp data of Krusche et al. p^. In addition, the fit included -kN 

''pi \ amplitudes of Arndt et al. |2J], since the irN and rjN channels are so strongly coupled. Using the notation for the 

' , ■ T-matrix: T~^ -\- iq^i — 1/a + ^Qn + sq^, - with g^ being the 77-momentum in the rjN center-of-mass - resulted in the 

O ' ryA^-effective range parameters 

:3 a(fm) = 0.75(4)-fi0.27(3), ro(fm) = -1.50(13)-i0.24(4) and s(fm3) = -0.10(2)-i0.01(l) 

^ A later paper involving the present authors |g^ developed this formalism to look for complex poles in the S- 

J> . matrix. In Ref. |£l| the only place, where the r]N amplitude T^^ was directly checked against the data, was under 

l^ ' the assumption that it was proportional to the photoproduction amplitude - namely - 



<^ ' _/.^„^ _ Arnu^-,\ /77^7T^r^7~r77I~^r^ „.-u _/.,„^ /_^.,„N ^7 



T{jr,)r = A{Phot)^J{Re T,„Y + {Im T,„Y , where t(777)^ = \h{Tl) 



Airq,, 



and A{Phot) was a free parameter adjusted alongwith other parameters to ensure a good fit to the combined data. 
However, in Ref. p5| an alternative assumption, that A{Phot) could also be energy independent, was also considered, 
namely, 

A{Phot) ex a(l -I- iT^^) + — T^„ (1) 

- a form that had proven successful in Ref. |2q| when describing pion photoproduction. This led to two distinct 
solutions called GWll and GW21. The former, with the energy independent A{Phot), resulted in parameter values 
for the complex poles that were close to those given by the Particle Data Group [E^ and are here referred to as 
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conventional (C). On the other hand the use of Eq.|l| yielded very unconventional (UC) values. In particular, the 
couphng of the ttN channels to the resonances at 1535(20)MeV and 1650(30)MeV generated widths of 354 and 133 
MeV respectively, whereas the GWll solution gave 152 and 150 in agreement with the PDG values of the resonance 
widths of 150(5) and 145 - 190 (all in MeV). Of course, it must be remembered that the PDG values are for poles in 
the T-matrix and not the X-matrix and that the nomenclature C versus UC is based on the underlying assumption 
that the T-matrix and iiT-matrix poles are closely related. Therefore, if, with an improved data basis, the GW21 
type of solution becomes the preferred solution ~ not the case at present - then it could indicate that the K-vaaXrvx. 
formulation generates strongly energy dependent widths. These differences in the widths are also reflected in the 
effective range parameters (a, tq, s). As expected, GWll gives values for a, tq, s that are very similar to those of 
the earlier fit in Ref. [£1|, whereas those from GW21 are very different^ being 

a(fm) = 0.21+10.30, ro(fm) = -2.61+16.67 and s(fm3) = -0.39-13.67, 

where the large values of the Imaginary parts give support to the possible need for energy dependent widths. This 
wild behaviour with such large imaginary components shows that the effective range expansion, in contrast to GWll, 
is poor for GW21. The results from Ref. ^^ and GWll are within the error bars quoted in ||2^]. However, they are not 
identical, since a slightly different selection of data and their associated error bars was used. In discussions concerning 
quasi-bound states, it is the value of T^^ in the unphysical region a few lO's of MeV below the r/A^-threshold that is 
relevant. In Figs. |l| a) and b) this is shown for both GWll and GW21. There it is seen that the GW21 amplitudes 
are not only smaller than those for GWll at the threshold but also they tend to drop faster over the range below the 
threshold that is important for quasi-bound states. The conclusion from (21| is that the 77 A^- effective range parameters 
seem to be crucially dependent on the form of A{phot) - a point that cannot be resolved within the restricted form 
of the if-matrix model presented there. 

In a third paper by the authors 0, the ii'-matrix model was extended so that, not only are the ttN and rjN 
channels treated explicitly - as in Ref. [£1[ -, but also the jN channel. This means that the photoproduction 
reactions jN -^ ttN and jN -^ rjN can be treated on the same footing as nN -^ ttN and ttN — » ijN. This avoids 
the need to choose a form for A{phot) to relate the 7iV -^ rjN and rjN -^ -qN processes. The penalty that is paid for 



this extension is that four more parameters are needed. In the two channel form of Ref. |21|, the basic parameters 
[£'o,ij 7Tr,77(0, 1), Bij] were those describing the various X-matrices - namely 



77r(0) , 7^(1) ^ _ _y7^(0)7^ , p, ^ Ir, 

p p n" B^ ' '^^1 ^ w w '^'' ' '''' ^ w w 






The £^0,1 are the positions of poles that in a "conventional" model should be near the energies of the 5*- wave ttN 
resonances A^(1535) and A^(1650). The 7,r,77(0, 1) are channel couphng parameters that are related to the widths of 
these resonances. Again these widths are thought to be more or less known when data is analysed by a "conventional" 
model. However, as mentioned above, less conventional models can lead to widths that are quite different. Finally the 
Bij are energy independent background terms and are purely phenomenological. In the extended iiT-matrix model of 
Ref. B the introduction of an explicit 7A^ channel requires the additional A'-matrices containing four more parameters 
7^(0, 1), Byr,, ByT,: 



VT^IOh^ , ^ ^ 77(0) , 77(1) r^ _ ^77 (0)7^0) , V77( 1)741) 

P P '^'' ' '^^ ~ TP fP P P ' ''^ ^ P P P P 



It should be added that the ttttN channel was also included implicitly in the formalism at the expense of two more 
parameters 73(0, 1) that are related to the tttt branching of the two resonances. These lead to four more AT-matrices: 



^ 73(0) 73(1) j^ _ ^7^0)73(0) ^7^1)73(1) ^ _ VtVT^ 

33 T-n rp rp TT* ' ^^ Tp Tp rp rp '^ ^3 rp rp "> 

J. _ ^77(0)73(0) ^77(1)73(1) ,. 

''^ Ea-E Ei-E ' ^ ' 

The ttttN channel is then removed by making, for example, the replacement 

^""Ao- A +^"" + '1-793^^33' ^^ 



^The authors wish to thank Bengt Karlsson for pointing this out. 



where q^ is a three body phase space element. It should be emphasised that this treatment of the ttttN channel is 
not an approximation. It is purely for convenience. More details can be found in Ref. f^. The outcome of that 
work was essentially two distinct solutions (A and B) depending on the number of 'jN -^ ttN data points taken 
from Ref. [^. Using only 32 points with Ecm < 1550 MeV gave a{fra) — 0.87 + iO.27 for solution A, whereas with 
48 data points and Ecm < 1650 MeV gave aj fm ) = 1.05 + «0.27 for solution B - values not very different to the 
a{fm) = 0.75(4) + i0.27(3) of our earlier work JIll. 

The two solutions given in Figs, n^ a) and b) may be understood in terms of a simpler model. There are three 
basic components of iiT^^ in Eq. ^. The pole term with Eq — E = /^ \s due to the A'^(1535) and 7,, is the appropriate 
coupling constant. The second term i?^^ represents background interactions due to an unknown potential in the 
elastic channel. The final term, which we now call iKine is a contribution from the inelastic (mostly ttN) channels. 
It contributes the absorptive part to if^j, and this is stressed by the factor i. Close to the threshold we can write 
A — Eq — Mri — Mn — Ekin = Aq — q^ /2mrfN, which induces some energy dependence and determines the effective 
range expansion. The scattering length now becomes 

Anrj{E) = ^n-i- + B,„ + iKine, (6) 

and the energy dependence of the K matrix is given by 

1 1 



which reduces to the usual effective range expression 

1 _ 1 g'7^ 



Kr,^[E) Ar,^[Q) Al„Al2m^N' 



(7) 



(8) 



provided A^^(O) is not small. Now, the solution of Fig. |l|a) corresponds to a dominance of the pole term over the 
background term -B^^. This induces a large scattering length and also a large effective range, so that Eq. ^ is valid. 
Nevertheless, the alternative expansion in Eq. also converges very fast in a broad energy range - see Fig. 2 in Ref. [£1| . 
On the other hand, the solution of Fig. ||b) corresponds to a strong cancellation of the pole and the background terms 
Jtj-^ — h Bjfjf w 0, so that Re Arjn ~ 0. In this case the effective range expansion in Eq. is not appropriate. A 
better parametrisation is obtained with simply an expansion of K^jniE) instead of l/Krjr]{E). Actually the zero in 
the scattering length occurs in the subthreshold region. This solution resembles to some extent the case of elastic ttN 
scattering. It should be added that those analyses that extract the rjN scattering length (a) from a many hadron 
reaction - such as pn ^ dr/ Q, Tyd — > ryd 0] or pd — > ^HeT] |]15[ - are, in fact, discussing a at an energy of 10 - 20 
MeV below the threshold. As can be seen from Fig. ^) this results a value of T^,, that is substantially smaller than 
the threshold value. Therefore, in those cases where an energy extrapolation is not performed, one should not be 
surprised, if the values of a quoted as being appropriate for such many hadron reactions are smaller than from those 
reactions involving at most two hadrons - as in the ii'-matrix description by the present authors. 

Since the extension of the AT-matrix model to three explicit channels removes the need to choose a form for A{phot), 
it is now possible to check whether or not the unconventional solution GW21 in Ref. |2^ is still acceptable. Now 
the most striking difference in the parameters in Table III of Ref. ||2^ is seen in the i?,,,,, which is 0.11 for GWll 
and -0.83 for GW21. There is a direct correspondence between a and B^^ - the more positive U^j, gives the more 
positive a. However, in spite of trying many different starting values in the Minuit minimization program - especially 
starts with large negative values for the background term B,y^ - the GW21-like solution could not be found. This 
suggests that the energy dependent form of A{Phot) in Eq. |l| is not a good assumption. However, this is not the 
case. The extended iiT-matrix model now supplies the input for checking both assumptions for A{Phot) - the energy 
independent form and the one in Eq. || - and in both cases the forms are found to be very nearly proportional to T-y^ 
- the differences being very small. This now suggests that the appearance and disappearance of the unconventional 
solution could be due to two almost degenerate minima in the parameter space - a feature that will rise again later. 

Since the publication of Ref. Ig] there has become available more experimental data: 

a) Firstly, the GRAAL collaboration has extended the 7p — > -qp data of Krusche et al. |g^ from the range 1486 < 
Ecm < 1539MeV upto about Ecm = 1700MeV [g^. This new data not only gives total cross sections but also 
differential cross sections, which do not concern us at this stage. When the new total cross sections, at the higher 
energies, are superimposed on Fig. 3 in Ref. g depicting the isT-matrix extrapolations beyond Ecm = 1539MeV, it is 
found that solution A corresponding to a{fm) = 0.87 + «0.27 is favored over solution B. 



b) Secondly, for the ttN -^ rjN reaction there is now new data |30| and also a reassessment by Strakovsky et al. |3l|] 
of the earlier review by Clajus and Nefkens p^]. The main difference compared with the earlier set is that, in the 
reduced cross section 



cr{TTf]),. = (j{TrT]) — , 

there is now a distinct structure over the pion laboratory momentum range 700 < p^ < 720 MeV/c - see Fig. g. In the 
earlier review [g2|, a{-K'q)r was almost constant with a small monotonic decrease upto about 770 MeV/c. If the data 
points in the range 700 < p^r < 740 MeV/c are ignored, then the Minuit fitting procedure always results in a unique 
solution independent of the parameter starting values. This appears to be a robust result. However, if the points in 
this range are included, then the resultant effective range parameters depend on the starting values - particularly on 
the background term B^^. Broadly speaking, if the starting value of B^^j is chosen to be positive, then the resulting 
solution always has Re a « Ifm and the resonance widths are conventional. However, for negative starting values of 
_B^^, the unconventional solution emerges with Re a « 0.2fm and resonance shapes very different to the Particle Data 
Group [g^. The resulting x^/dof are very similar for the two cases being respectively 1.036 and 1.029 with 155 data 
points and 12 parameters. This clearly shows that, within the present ii'-matrix model, there are two minima in the 
parameter space. However, as seen in Fig. g, the reason why the unconventional solution emerges, as the data points 
in the range 700 < Ptt < 740 MeV/c are incorporated into the data base, is that the featureless curve for Re a k, 0.2fm 
is lower than the equally featureless curve for the higher values of Re a. Even so, neither solution is able to account 
for the structure now seen in the reduced cross section. This result can be interpreted in several ways: 

a) In the iiT-matrix model, some effect is omitted in t:N — » rjN channel in the present form of the parametrisation. 
But it is difficult to image anything that could have such a dramatic effect so near the threshold. 

b) The results in Refs. JS^,^^ need to be checked, since if this structure is indeed confirmed, it indicates the presence 
of some hitherto unexpected mechanism not seen in the other channels. 

In Table II, the parameters corresponding to the conventional (C) and unconventional (UC) solutions are compared 
with each other and the values given by the Particle Data Group p^. There it is seen that the two solutions are 
qualitatively the same as GWll and GW21 in our earlier work ES]. Again the distinguishing features between the 
two fits are the values of the energy of the X-matrix pole and the corresponding widths [Eq and r(1535) in Table 
II]. Of course, this is not reflected in the resonance position and width defined as the Real and Imaginary resonance 
energies from the T-matrix. In fact, they are almost the same for GWll and GW21 with the values Ep =1514 versus 
1509MeV and r/2 =90 versus 82MeV respectively. The difference between the two solutions is the non-Lorentian 
shape of the resonant T-matrix in the GW21 case. This is seen clearly in eq.(18) of Ref. [||. and follows the same 
analysis that leads here to Eqs. um 

Even though the overall x^/dof are virtually the same for the two solutions, the individual x^ /dp for each type of 
amplitude are much different - as can be seen in Table III. 

One of the authors (S.W.) wishes to acknowledge the hospitality of the Helsinki Institute of Physics, where part of 
this work was carried out. The authors also thank Drs.T. Morrison, D. Rebreyend, F. Renard and I. Strakovsky for 
useful correspondence concerning their unpublished data.. 
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Figure Captions 



Fig 1: The r/N -^ r]N amplitude T{j]N) — Tjjri in fm as a function of the center-of-mass energy Ec.m. in MeV. 
Solid line for Re T and Dotted line for Im T. a) is for the conventional(C) solution GWll and b) is for the 
unconventional(UC) solution GW21. 

Fig 2: The reduced nN -^ rjN cross section a{7Tr])r in mb as a function of the Tr-momentum in MeV/c - see Refs. 
pO| and |pi[. Sohd hne - conventional solution(C), Dotted line - unconventional(UC) solution. 



TABLE I. A selection of r^N- scattering lengths a appearing in the literature. 



Reaction or Method 



Scattering Length(fni) 



Isobar model 
Isobar model 

Isobar model [ HI 
pd -^ ^Her] M 
Coupled T-matrices 
Effective Lagrangian 
Coupled _R'-matrices 
J7d -> J7d § 
Coupled ii'-matrices 
Coupled T-matrices 
Coupled T-matrices [ Lij 
Coupled ii'-matrices 




0.25-fi0.16 
0.27+10.22 

<0.3 
0.51-M0.21 
0.55(20)-Hi0.30 
0.621(40)+i0.306(34) 
0.68+10.24 
0.75(4)+i0.27(3) 
>0.75 
0.87+10.27 
0.886+10.274 
0.98+10.37 
1.05+10.27 



TABLE II. The optimised parameters defining the ii'-matrices for: 
tional(UC) solution compared with the Particle Data Group pM. 



a) a conventional(C) solution and b) an unconven- 





C 


UC 


PDG 




C 


UC 


PDG 


-Br, 77 


0.4225 


-0.4769 




r(1535, Total){MeY) 


125.5 


267.0 


100-250 


S,rr7 


-0.0376 


-0.0282 




r;(1535, br) 


0.590 


0.578 


0.30-0.55 


£;o(MeV) 


1529.7 


1580.5 


1535(20) 


7r(1535, br) 


0.346 


0.358 


0.35-0.55 


£;i(MeV) 


1681.3 


1681.0 


1650(30) 


r(1650, Total){MeY) 


165.6 


164.6 


145- 190 


-D717 


0.00223 


0.00707 




7r(1650, br) 


0.732 


0.731 


0.55 - 0.90 


Bjn 


0.00250 


0.00484 












77 


0.000166 


0.000240 













TABLE III. Comparison of the X'^/dp for the separate amplitudes. In the tvN 
refer to the contributions from the Real and Imaginary terms respectively. 



nN and 7A'" -^ nN cases, the two entries 



Solution 

C 

UC 



ttN -> riN 
0.981 
1.017 



ttN -^ nN 
1.281, 0.803 
1.206, 0.808 



'yN ^rjN 
0.869 
0.893 



-yN -^ -kN 
0.556, 1.403 
0.538, 1.334 



Total 
1.036 
1.029 



FIG. 1. The r]N — > r]N amplitude T{riN) = T^^ in fm as a function of the center-of-mass energy Eq.m. in MeV. Sohd line 
for Re T and Dotted line for Im T. a) is for the conventional(C) solution GWll and b) is for the unconventional(UC) solution 
GW21 
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FIG. 2. The reduced ttN — > rjN cross section aijvri)^ in mb as a function of the 7r-momentum in MeV/c 
[pl[. Sohd line - conventional solution(C), Dotted line - unconventional(UC) solution. 
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